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Phonon polarization in a magnetic field is analyzed in isotropic model. It is shown, that at presence 
of spin-phonon interaction phonon possess circular polari-zation which causes the appearance of heat 
flux component perpendicular both to temperature gradient and magnetic fleld. 
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Recently a novel phenomenon was found. It 
is the detection in dielectric of a heat flow perpendicular 
both to magnetic field B and temperature gradient VT - 
so called Phonon Hall Effect (PHE). Such fiow is caused 
by spin-phonon interaction (SPI) of phonons and para- 
magnetic ions r&+^. The theory of PHE was considered 
in works j3| and jl] (the critic of work [3| see in 

However not only SPI leads to magnetic field depen- 
dence of dielectric heat conductivity. As it is known, in 
molecular gases magnetic field affects the transport phe- 
nomena because of the molecular rotary moments M, 
precessed in a magnetic field, dependence of coUision 
cross-section (Senftleben-Beenakker Effect Q). 

We expect, that the phenomenon similar to PHE in 
Tb3Ga50i2 can be found out in molecular crystals where 
is a component, consisting of the molecules with rotary 
degrees of freedom (RDF). Typical representatives of 
such substances are criocrystals To magnetic field de- 
pendence of heat conductivity can lead a RDF-dependent 
phonon scattering on molecules. The suitable theory is 
easy to construct via generalizing Senftleben-Beenakker 
theory on a case of heat transfer by phonons. In approach 
where RDF are considered quasi-classically such theory is 
considered in Q • In this paper we have another approach 
to this problem, and we show that in crystals with RDF 
there is a mechanism close to SPI in ionic crystals which 
leads to PHE through renormalization of acoustic waves 
because of lattice vibrations interaction with molecular 
rotation. From symmetry considerations describing such 
interaction Hamiltonian can be written down in similar 
to SPI form: 



Hi 



Pn]) . 



(1) 



Here M„ is the rotary moment of the molecule in a cell 
n, [u„ X p„] - total orbital moment of oscillating atoms 
in a cell. Speaking generally, molecular crystals have 
a complex structure containing of many particles in an 
elementary cell. But at low temperatures heat is trans- 
ferred by long- wave elastic waves, when all particles of 
a cell oscillate with the same amplitude u„ and velocity 
v„. Therefore instead of a multicomponent crystal it is 
possible to consider a crystal with one atom in a cell with 
mass TO, equal to total mass of particles in a cell and one 
rotary moment. The value of g in case of SPI was es- 
timated in numerous works on spin-lattice relaxation j3| 



and [3| . In case of acoustic vibrations interaction with ro- 
tary degrees of freedom factor g we shall consider as small 
phenomenological parameter having the same order as 
for SPI when the molecule has uncompensated electronic 
moment, as in 02- We shall use system of units where 
ks = I, h ^ I. For short, interaction ^ shall be also 
called SPI in the case of crystals with rotary degrees of 
freedom. Almost identical form of Hamiltonian describ- 
ing phonon interaction with internal degrees of freedom 
in ionic and molecular crystals (components of opera- 
tor M are replaced with Pauli matrices in ionic crystals 
with quasi-doublet structure of the bottom levels) leads 
to qualitative similarity of transverse thermal conductiv- 
ity theories in both cases. But from the physical point of 
view interaction mecha-nisms are far from each other. 

Usual dielectric heat conductivity weakly depends on 
detailed properties of acoustic vibrations. As it is shown 
in PHE is caused by occurrence of phonon elliptic 
polarization at presence of a magnetic field and SPI. But 
character of acoustic branches polarization essentially de- 
pends on symmetry properties of a dynamic matrix. Un- 
like the general case considered in \^ , where all three 
vibration modes are nondegenerate almost at all direc- 
tions of wave vector k, in present paper the model of 
vibrations in isotropic medium where transverse acoustic 
modes are degenerate will be considered at all k. Though 
this model is reahzed only when atoms with a large in- 
teraction radius interact, in this model dynamic matrix 
has rather simple form. In isotropic medium model we 
shall find a heat fiow in a [B x VT] direction as for crys- 
tals with rotary degrees of freedom, and for ionic crystals 
with quasi-doublet structure. The first result is interest- 
ing as a prediction of PHE presence in molecular crystals, 
and the second is correction of work j^l- 

In an external magnetic field the average value of the 
rotary moment (M) = (M„) is nonzero, and crystal elas- 
tic vibrations are determined by Hamiltonian renormal- 
ized taking into account Jl]): 
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where 



hn = ^pI -\Y. ^Dt,«'-9 (M) [u„ X p„] . (3) 

n' 

It is essential, that the effective couphng constant g (M) 
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is proportional to magnetization of particles and it dis- 
appears when someone switches off the magnetic field. 

Using ([3]) it is possible to write the quantum energy 
continuity equation and to deduce the density of a heat 
fiow formula (|3, d, d) 



fH-^Y.''Dt'«'- (4) 
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We emphasize, that in this expression term v'^ is a veloc- 
ity of this ion, but not an impulse of an ion divided 
by mass 

< = dtK = = Pn/m - eabcg (M") <. (5) 

Let's turn in (0]) to impulse representation using 
J2 r'^Df exp (zfcr) = iVlOf. 



We find 



^H=^^T. i^lDf)uLv\,,. (6) 



(s is the mode number) Let's introduce decomposition 
of displacement vectors and velocity vectors for normal 
vibration modes 



(7) 



It is very important, that usual form has a veloc- 
ity decomposition instead of an impulse Pi — mvi — 
gm (ui X M) as authors of Q supposed. We substitute 
u and D in ([6]) for ([7]) and average by state, which is di- 
agonal on phonon numbers. Then we discard anomalous 
averages {aksO,-ks'} and (alj,g,a^^) and we come to heat 
fiow in a [-B X VT] direction: 
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In isotropic body PHE is characterized by one factor (as 
a Hall constant). Without loss of generality we choose 
magnetization on Z-axis, a temperature gradient - on X, 
and a heat fiow - on Y . Using this expression in Q a 
PHE amount for acoustic vibrations of a general view has 
been calculated. In present paper we consider an elastic 
vibrations model in isotropic body where the dynamic 
matrix in zeroth-order approximation D'^ has a form 

Df = cl6''''e +wkah, (9) 

with degenerate branches = loq = c^k^ and longi- 

tudinal branch = Wq -I- A, A = wk'^. Here the third 



independent term like (5"^fc^, which in a cubic crystal re- 
moves degeneration is discarded. 

For isotropic body model considering orthogonality of 
polarization vectors e^*e^^, = 5ss' we have 



+ M(h^Z}4.' + 4: (*:kL.))1 «.<■/..'>- (10) 

We see that ifTOj) do not contain factor g evi- 

dently. The heat fiow magnitude changes only because 
of renormalization of phonon spectrum and polarization. 
Everyone can verify that diagonal terms in lfTO|) do not 
lead to heat fiow in [B x VT] direction. PHE is described 
by off-diagonal terms in lfTO|): 



(11) 

Particularly, transverse modes correlated motion gives 



+ C(kefe-)(4+«fc-)}- (12) 



Expression (fT2l) essentially depends on phonon polar- 
ization vectors. First of all it is necessary to establish 
their form in consideration of SPI. Hamiltonian ^ leads 
to the dispersion equation 

^ks^ks = Df'e^s, (13) 

where 



(14) 



[G) ^ [Df i-G)) , {bI,bus') = S,s'. (15) 

SPI contribution to dynamic matrix is an imaginary 
antisymmetric tensor. In model ^ dispersion equation 
takes a form 



W^e° = Ar(fce) + i[ex G]" 



(16) 



Here W =^ uj"^ - ujI, \ = wk^, k = k/k = 
(sin6'cos<y5, sin0sin(y9, COS0) and G = (0,0, G), (^Gfc^ — 
Gcosd, Q = Gsin6'. Let's use orthonormal basis vectors 



l)k, 2)m = (G - [Gkj k)/Q = 
= (— cos 9 cos (f, — cos 6 sin cp, sin 6) , 
3)n = k X 771 = Q^^ kxG = (sin — cos (/s, 0), 
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From (fTBl) we find vector e components 



Q 



and eigenvalue equation 
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(17) 
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Dispersion equation (fT3|) determines polarization vec- 
tors accurate within phase. We set phase of parameter 
^ equal to zero. Parameter absolute value is set by nor- 
malization: 



(19) 



At inversion (k,m,n —k^m^—n) polarization vector 
(fTTl changes its sign as a polar vector. Complexity of 
vector ifTTl means that it is eUiptically polarized. 

Expressions ifTT)) . (fTSj) . ifTO]) are correct with any G. 
Then we assume phonon and internal degrees of freedom 
of ions (molecules) interaction to be weak G <C A. In this 
case we get one longitudinal mode W\\ = A, e = fc from 

For transverse modes in zeroth-order approximation W 
is equal to zero, and in consideration of G is equal to 



W„ 



1 



and describes their splitting 



(20) 



1-J (Qy2Xf + (GkY . (21) 



Splitting is minimal on equator (|cos0| < G/A, A 
Q^/ (2Aajo))- For other k directions we have 



Wrj^rj Gk 



G |cos 6*1 /luq. 



(22) 



Real and imaginary parts of each e,, in zeroth-order 
approximation mutually perpendicular end equal to each 
other. 
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[rj (sign cos 9)rh — in] 



(23) 



In linear approximation we find 



We also need following expressions: 
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(sin sin (/s) 



4AW,, 



(24) 



(25) 
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+ icos(p\. (26) 



Expression (fT2l) contains projection of polarization vec- 
tor onto wave vector direction ^e^A:^ which is Hnear on 

SPI. So, off-diagonal density matrix {a^gaks') is quite 
enough to calculate in zeroth-order approximation on 
SPI. 

As is known, the nonequilibrium part of density matrix 
(diagonal on modes) in tau-approximation is equal 



1 



/p = {a+ap) -Np = -— (CpV) 

i Lpp 

1 



Np (1 + Np) (ijpCp) VT. (27) 



Here and afterwards, we use p = ks for short, l/ilpp is a 
p-phonon relaxation time. Furthermore, in obvious cases 
we shall leave out an index k. 

The corresponding energy flow is parallel to VT and 
leads to heat conductivity coefficient 

>f^^ ~ {cny^ , (28) 

where c, - average values of Cp, ilpp. 

Off-diagonal part of density matrix can be expressed 
by fp constructing generalized Boltzmann equation. 
Phonon scattering has many channels: anharmonicity, 
resonant scattering with molecular multiplet structure 
excitation, scattering with rotary degrees of freedom con- 
tribution, impurity scattering. In all cases the answer has 
similar structure. In Q it is shown, that 
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(29) 



This means transverse phonons have circular polarization ^^ere Jp„ is a Hermitian matrix which differs mainly 



Polarization vector (|23|) direction changes stepwise 
crossing equator. However, if we change transverse mode 
numeration and write W,^ = -qGk instead of l(22|) then 
projection of (fT7|) on the m-axis will be constant. It al- 
ways happens when crossing levels. The magnitude of 
(fT2l) does not depend on numeration. 



from usual collision integral in that there is product of 
the phonon scattering amplitudes instead of a square of 
the scattering amplitude module (Fermi's golden rule). 
In tau-approximation this expression has a form of 



) 



2 {Ujp - LUq) 



(30) 
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find accurate within numerical coefficient 



^,._iMV|;«,(l + A.,).-iMI ,38, 
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Dividing it on longitudinal heat conductivity coefficient 
~ r3 we find a Hall angle 



where VL*^ {up) = (ujp) and (a+aq)* = (a^ap) are Consequently, accepting w ~ and 
Hermitian. 

The effective relaxation frequencies form depends on 
scattering mecha-nism. In [4] corresponding expressions 
for resonant scattering and anharmonicity are shown. In 
case of potential impurity scattering 

^qp (wp) = 27r^^ Aa^gp^ i^P " ^g) ■ (31) 

g 

It is essential, that the transition amplitude Apq — 
ep*Ap^eq is proportional to polarization vectors. It guar- 
antees the absence of a stream (fT2|) dependence on gen- 
eral phases of phonon modes (phases of parameter ^ in 
(fT?!) '). The Aqg polarization-dependence is automatically 
transferred on relaxation frequencies ^Iqp = eq*ftqpep, 
where tensor = (^pg)* does not depend on exter- 
nal polarization vectors. In isotropic body model we can 
neglect mode splitting for transverse modes in numerator 
of Il30l) 
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The second channel of Hall heat conductivity is caused 
by occurrence in temperature gradient presence a corre- 
lated transverse and longitudinal phonon motion which is 
characterized by off-diagonal density matrix component 
131: 



(a+aa) 
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[^3v M fl + ^7)3 (^3) /s] 



(39) 



(32) 



Here 1 = (k, _L) , 3 = (k, ||). This density matrix (and its 
Hermitian conjugation) forms Hall heat ffow (see (fTTj) ) 



When one of modes is longitudinal from ((3Ql) we have 



{a^a,j) = ^[^3»? (ws) /s + f^r,3 (^0) fl 



(33) 
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where e = ui^ — ujo = \/^q + ^ "-^o . 

Let's use the formulas above for calculation of tranS' 
verse heat conductivity coefficient ^c^^. Using ((32|) in lfT2] 
we find 



+ ef (fcei) (a+ai)] (40) 



W \ - 

2V ^ 

k 



Rc{-e^l (kefc+) in*_^ (cjq) ^ + 

+ e^;(ke,_)zf7_+(c.o)^}. (34) 

In (|27|) function f± is proportional to ~ cos(p, and 
nonzero result is given only by imaginary part of (|26l) : 
e^* = (...) - iS^^cosip. From ^ and dSl]) we have 

/■V\ W yr-^ fj_kQ 

xRe{-(l + ^)n%+(l + ^)n^+}. (35) 



The matrix (|33| is proportional to ^7\ ) and cosLp. Be- 
sides, e^* = —i£,ri cosLp and e|* = sin sin 93. So, only the 
first term in square brackets gives a contribution to inte- 
gral. In this case ^^63^ ~ 1 and integration element in 
(|40l) does not contain series expansion parameter of SPI. 
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— ) cosip^ifi Refi3i + 
ivi 2e 

/3Rcrii3(a;3)) (41) 



Using (|27|) and integrating we estimate longitudinal mode 
contribution to Hall heat conductivity 
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r2 Re 1^31 



Co 



(42) 



We see, that only second-order terms in (|24|) give the 
contribution to PHE. Using concrete form of parameters 
A, A, Q, W^, G we find 



" (c|| - Co) ripp 

In this case the Hall angle has an order of magnitude of 
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T 



(43) 
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Then we integrate using ([23 

g (M) w ^ wgc 



x" = — - 



V 



k 



2 Re 17-+ 

(sm6'cos<y9) 

xNp{l + Np). (37) 



In conclusion we compare the estimations of heat con- 
ductivity coefficients caused by transverse and longitudi- 
nal mode correlation k^^^ (|42l) , transverse modes (|38l) , 
and estimation x^j^j for ion crystal with nondegenerate 
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(44) 
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We see, that estimations for transverse modes in isotropic 
model and cubic crystal coincide accurate within average 
polarization caused by rotary moment substitu-tion for 
pseudospin polarization. It is essential, that in both cases 
frequency de-pendence is absent as it should be for heat 
flow component in [B x VT]. direction. The longitudinal 
mode l(42|) contribution has qualitative other form. It 
does not depend not only on frequency but on coupling 
constant too, as it is in common Hall Effect. At high 
temperatures (T ^ g) channel plays a main role. 
If there were crystals where even approximately there 
were acoustic vibrations with ellipticity close to circular, 
then PHE experimental observation would become rather 
simple problem. 

The flrst theoretical work devoted to PHE 01 contained 
a mistake at using relation ^ between particle velocity 
and their impulse at presence of SPI. Nevertheless it is 
interesting to compare a estimation in (3] with our 
formulas l(42|) . l(38|) . as in [3] calculations were performed 
in isotropic model like in present paper. If to not pay at- 
tention to numerical coefflcients which are obvious excess 
of accuracy, they re-ceived a result 



which coincide with fSS]) and main result of This con- 
currence is caused not by of our theory equivalence, and 
only trivial consequence of similar dimension of results 
and that all of them are received in linear approach on g. 
So, it is necessary to emphasize, that the basic result of 
the this work (and Q) is not the estimation of , and 



the establishment of that fact, that PHE is caused by 
joint action of two equally important factors - 1) phonon 
elliptic polarization caused by SPI, and 2) induced by 
a temperature gradient the correlated movement of two 
phonon modes with off-diagonal density matrix forma- 
tion. 
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